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ABSTRACT 

We present results on the evolution of the intrinsic scatter of black hole masses consid- 
ering different implementations of a model in which black holes only grow via mergers. 
We demonstrate how merger driven growth affects the correlations between black hole 
mass and host bulge mass. The simple case of an initially log-normal distributed scat- 
ter in black hole and bulge masses combined with random merging within the galaxy 
population results in a decreasing scatter with merging generation/number as pre- 
dicted by the Central-limit theorem. In general we find that the decrease in scatter a 
is well approximated by crmerg(w) ~ (7i„i x {m -\- 1)"°/^ with a = 0.42 for a range of 
mean number of mergers m < 50. For a large mean number of mergers (m > 100) we 
find a convergence to a = 0.61. This is valid for a wide range of different initial distri- 
butions, refill-scenarios or merger mass-ratios. Growth scenarios based on halo merger 
trees of a (100 Mpc)^ dark matter ACDM-simulation show a similar behaviour with 
a scatter decrease of a = 0.30 with typical number of mergers ttt, < 50 consistent with 
random merging (best matching model: a = 0.34). Assuming a present day scatter 
of 0.3 dex in black hole mass and a mean number of mergers not exceeding m = 50 
our results imply a scatter of 0.6 dex at z = 3 and thus a possible scenario in which 
overmassive (and undermassive) black holes at high redshift are a consequence of a 
larger intrinsic scatter in black hole mass. A simple toy model connecting the growth 
of black holes to the growth of ACDM dark matter halos via mergers, neglecting any 
contribution from accretion, yields a consistent Af, — Msuige relation at z = - if we 
assume the correct initial relation. 
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1 INTRODUCTION 

There is growing observational evidence that most if not 
all bulges and elliptical galaxies hos t a supermassive black 
hole (SMBH) at the p resent time (|Magorrian et al.1 [l99i . 
iGenzel fc Eckartl ll999l. Furthermore there exists a strong 
correlation between black hole masses and properties of the 
host galaxy, e.g. the host bulge mass and the bulge velocity 
dispersion JH aring fc Rix' '2004', 'Ferrares e fc MerrittI |2000| . 
iGebhardt et al . 2000, Tromaino ot al. 200^) a nd possibly the 
host h alo ( Ferraresell2002l' ) in nearby galaxies. iHaring fc Rixl 
(1200J) find the relation between black hole mass M. and 
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bulge mass Meuigc to be: 

log(M./M0) = 8.20 + 1.12 X log(A/Buigc/lO"M0). (1) 

The correlation between black hole mass M. and velocity 
dispersion at can be written as 



log(Af./M0) = alog(cr./200kms"^) + b, 



(2) 



where a is the slope and b is the zero point. In the literature 
the values fo r the slope vary betw e en a = 3.68 and a = 4.86. 
For example iFerrarese fc Merritti (|2000l) give values for the 
zero point b = —2.9 and the slope a = 4.80. Ge bhardt et al 



2000) claim a smaller slope of a = 3.75, T remaine et al 
2002) find a = 4.02 and more recen tly IFerrarese fc Ford 



(2005) estimated a slope of 4.86 and iGrahamI l|2008l ) one 
of 3.68 for barless galaxies. Concerning the intrinsic scat 
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Figure 1. Comparison of the intrinsic scat ter from observations 
llTremairie et alj|200l I CiJltekin et al.ll2009h and theoretical pre- 
dictions I IMalbon et al.l2007l ). From the 10-90 percentile spread in 
the model we have calculated the corresponding theoretical scat- 
ter. The shaded areas illustrate the error on the observationally 
estimated intrinsic scatter. 



ter a in black- hole mass, logM, , m ost studies agree that 
the scatter is not larger than 0.3 dex ([Gebhardt et al.l l200d : 
iTremaine et all l2002l : iNovak et all |2006^). Note, that with 
a* we describe the central velocity dispersion of a galaxy, 
whereas we characterize the in trins ic scatter as a. I n con - 
trast to iGebhardt et al.l (l2000l') and IXremaine et all l|2002h . 
a recent study bv lGiiltekin et al.l (|2009l ') obtains for Eq. [2]a 
slope a = 4.24 and an intrinsic scatter in black hole mass 
of a = 0.44 dex for a sample of 49 M.-measurements. For 
a subsample of early-type galaxies they find a smaller slope 
as well as a smaller scatter [a = 3.96, a = 0.31 dex) than 
for the full sample. For non-elliptical galaxies the intrin- 
sic scatter is larger with a — 0.53 dex. In the following 
the scatter a will always be given in dex. The existence of 
these tight correlations strongly suggest a co-evolution of 
the black hole and the bulge component of the host galaxy. 
However, the origin of these relat ions is uncertain and 
a subject of current research (e.g. IVolonteri fc NataraianI 
20091: |Pendl200i1: iBurkert fc Silkll200ll : ISpringel et al.ll2005l : 
Johansson et al.l 120091 ). Apparently, the origin of these re- 



lations can be connected t o the gas dynam ics in ma.jor 
galaxy mergers llMihos fc He rnouist 1996, Naa b et al.ll2006l . 
[Robertson et al.ll2006l . lHopkins et al.ll2008t ). In this scenario, 
the black holes grow significantly in gas rich mergers of 
disk galaxies and the remnants appear on the observed scal- 
ing relations. Subsequent, possibly gas poor, major and mi- 
nor m erging (Khochfar & Burkcrt 2003, Naab ot al. 2006, 
iKhochfar fc S ilk 2009, Naab ct al. 20 09|) conserves the rela- 
tion lISesanTet al... 2004. .Robertson et al.ll20o"6l. |Penell2007l. 



Hopkins et al. _ 2007al. Hopkins et al.l l2007d . ISpringel et al.l 
20051 . 1 Johansson et al.l l2009V 



Observationally the black hole mass relations are well 
constrained only in the nearby universe and it is unclear 
if and how they evolve with cosmic time. Several au- 
thors have found evidence that galaxies at higher redshift 
have a higher black hole to bulge mass ratio M,/Af Buigc 



than ellipt icals tod av IIMcLure et al. 120061: Treu et al.l 



Woo et al. 2008: Walter et al 



2004: Schramm et al 



2007; 



200S 



Peng et al. 2006; Greene et al. 2009; 'Nataraian fc Treisteij 
2009 ; Salviandcr ct al. 2007 ; Shields et al . 2006). For a sam- 
ple of Seyfert galaxies at moderate redshifts z < 0.1 the 
black holes are more massive by AlogM, ~ 0. 5 dex com- 
pared to the local blac k hole-bulge mass relation IjTreu et al.l 
l200'7l . lWoo erai]|2008l ). rSalviander et all (|2007l ) find at red- 
shift z ~ 1 an evolution of the M.-(j«-relatio n by 0.2 dex in 
black hole mass. At higher redshifts of z ~ 2 iMcLure et al.l 
(|2006t ) obse rve black ho l es 8 t imes more massive than ex- 
pected and iPeng et aU (|2006l 'l show that for 2 > 2 the 
M./AfBuigG-ratio is 3 — 6 times larger t han today. This has 

1120091) based on a lensed 

20081 ) find evidence for an 

10. At redshifts 



been confirmed bv I Greene et al. 



quasar sample. ISchramm et al.l 
excess in M, /A/buIkc at z ~ 3 of a factor of ' 
4 < z < 6 Shields ct al. (2006) obtain for black hole masses 
in the range of 8 < log(M./M0) < 10 a deviation from 
the present-day M .-Mbuige-relation of Alog(M.) ~ 2 dex. 
IWalter et al.l (|200J) report an even higher redshift object, a 
quasar at 2: = 6, whose black hole is about 20 times more 
massive tha n expected. Con s iderin g a present-day scatter 
ac cording to Trernaine et al.l (|2002l ) (a = 0.3) or according 
to iGiiltekin et all (|2009l ) la = 0.31 for ellipticals) all ob- 
served black holes at 2 J5 2 are outside the 2 — ct range 
of the present-day scatt er. Furthermore, rece nt observations 
( Alexander et al. 2008 ; iShapiro et al.l [20091) show also the 
existence of under-ma ssive black holes at high redshifts. 
[Alexander et al.l ([2008[ ) find black hole masses, which are 3 
times smaller than those found in comp arable massive galax- 
ies in the local Universe. The results in Shapiro ct al. ( 2009|) 
show that black hole masses at z = 2 are an order of mag- 
nitude lower than those predicted by local scaling relations. 
The most obvious explanation for the over-massive 
black holes at high redshifts are possible selection effects. It 
is more likely to detect the most luminous and most massive 
black holes at high redshift than less luminous ones. How- 
ever, the probability for finding a massive black hole in the 
mass range 10^ — 10^" A/© in the observed volume at z — 3 
(e.g. [Schramm et al. 20081 ) is extremely low a s estimates of 
the l ocal SMBH mass function from SDSS ([Benson et al.[ 
[20071 ) would predict no high mass black holes to be found in 
a similar volume even assuming no evolution in the SMBH 
mass function. Cosmic variance is very unlikely to be an ex- 
planat ion for the observed, massive black holes. [Lauer et al.[ 
((20071) point out that there is an additional bias which is due 
to different selection effects for high-redshift (e.g. black holes 
in high-z galaxies selected by nuclear activity) and local sam- 
ples (e.g. black holes in local galaxies selected by luminosity 
or velocity dispersion). They deduce that because of this 
bias M, will typically appear to be too large in a distant 
sample for a given luminosity or velocity dispersion. Some 
authors (e.g. Croton 2006) explain observed, over-massive 
black holes at high redshift with a shifted relation towards 
higher black hole masses for a gi ven bulge mass. He us es the 
Millennium ACDM-simulation ([Springel et al.[ [2OOII ) cou- 
pled with a model of galaxy formation, where galaxy merg- 
ers are the primary drivers for black hole and galaxy growth 
and explore an additional growth channel through which 
only bulges gain mass, e.g. the disruption of stellar galactic 
disks in major disk mergers. He argues, that if the bulge 
growth rate from such disrupted disks is not constant with 
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time, an e volution in the M . -Mbui ge-relation can occur. Fur- 
thermore, [Robertson et al.1 (|200q 'l find from simulations of 
galaxy mergers, that the M,-(j, -relation shows a slight red- 
shift evolution towards higher black hole masses for a given 
a, at higher redshifts, but the y predict no ev o lution for the 
M.-Meuige-relation. However, iHopkins et alj l|2007bl ) using 
again simulations of major galaxy mergers show that high 
redshift black holes will be more massive at a fixed bulge 
mass than expected from the present-day relation. They 
find an evolution towards lower black hole to bulge mass 
ratios with cosmic time which is driven by the fact that 
disks (merger progenitors) have characteristically larger gas 
fraction s at high redshifts. We want to point out that these 
studie s (|Crotonll200"^ : [Robertson et al.ll2006l : iHopkins et all 
l2007bl ) are consistent with each other, they are only focusing 
on different aspects of the evolution of the relation. However, 
these methods do not provide a sufficient explanation for the 
observed under-massive black hole masses at high redshifts, 
since these methods find to have black holes mainly lying 
above the median Af.-A/Buigc-relation at high redshift evolv- 
ing towards the relation with decreasing redshift. However, 
they also predict that the black holes go through a rapid 
growth phase before they end up above the relation; these 
objects might be under-massive. Furthermore, if a black hole 
at high z is already above the median M,-MBuigc-relation it 
will not end below the median relation at z=0 if only assum- 
ing merging and not gas accretion. Therefore, an alternative 
explanation for the observed high and low M,/Mbuigc-ratio 
at high redshifts could be the existence of a larger intrin- 
sic scatter in black hole mass, even assuming no evolution 
of the mean re lation with c osmic time, what would be in 
agreement with lLauer et al.i (2007 V 

In this paper we address the question: how does the 
intrinsic scatter in black hole mass evolve and change with 
time assuming that black holes grow only via mergers? An 
answer to this question is also important with respect to 
similarities and differences between the observed scatter of 
bla ck hole masses and predictions from theoretical mod- 
els. |Malbon_et_alJ (|2007|), using semi-analytic modelling, find 
that the present day scatter in black hole mass decreases sig- 
nificantly with increasing black hole mass. T his is in contrast 
to obs e rvations using th e full s ample s of e.g.lGebhardt et al.l 
l|2000l ). iTremaine et all (|2002l ). and iGiiltekin et al.1 l|2009^ . 
Here, the scatter appears to be independent of black hole 
mass. In particular at the high mass end the observed scat- 
ter is mu ch larger than the mo del predictions (see Fig. [T]). 
However, iGultekin et al.l (J2009|) demonstrate that the scat- 
ter for non-elliptical galaxies (typically at lower masses) is 
larger than for elliptical galaxies. This is, at least qualita- 
tively, in agreement with the predictions from lMalbon et al.l 
((20071). However, in particular at the high mass end the ob- 
served scatter is significantly larger than the model predic- 
tions (Fig. [1|. 

So far the time evolution of the sc atter in blac k hole 
mass has not been investigated in detail. iPend (|2007r ) deals 
with the evolution of the scatter, but he mainly focuses on 
the aspect of how the present day M,-Mbuigc-relation can 
form in a simple model applying random merging of galax- 
ies. He claims that the relation develops even if black holes 
and bulges are uncorrelated or incorrectly correlated in the 
beginning. This behaviour is supposed to result from an ini- 
tially exponentially decreasing SMBH mass function where 



minor mergers drive the objects towards the observed corre- 
lation. Furthermore, the scatter in black hole mass decreases 
with increasing merger number, and according to his results 
the decrease in scatter is dominated by major mergers. How- 
ever, a quantitative study of the scatter evolution was not 
presented, which is the main subject of this pap er. A f urther 
major difference between our work and that of lPeneJ ([20071) 
is that we include - besides Monte-Carlo generated random 
merging scenarios - merging as it is found for dark matter 
haloes in large scale cosmological simulations. 

In the following we will investigate the evolution of the 
intrinsic scatter assuming: 

• Simple random merging (section |2.2|) 

• Modified random merging (section 12.312.3.31) 

• Merging in ACDM-simulations (section [3]) 

We point out that random merging does not describe a full 
physical evolution process according to currently favoured 
structure formation models. In principle the model follows 
dry merging of galaxies and therefore is limited to high mass 
galaxies, since merging at the high mass end is assumed to 
be almost dry, so that gas physics and star formation do not 
play an important role and can be neglected in the growth 
processes. However, an advantage of using a simple model 
such as random merging is that we can study separate ef- 
fects on the scatter evolution, e.g. the infiuence of the initial 
mass distribution, of the merger mass-ratio or different refill- 
scenarios. Then we compare these results to merging accord- 
ing to currently favored structure formation models based 
on dark matter ACDM-simulations. Since we know that a 
significant contribution to black hole growth is caused by 
accretion, we will discuss this issue in section |4l 



2 MODELS FOR RANDOM MERGING 

2.1 Initial conditions 

In the following we describe two different initial distribu- 
tions of bulges including black holes as a starting point for 
random merging: a log-normal distribution and a Schechter 
distribution of bulges. 

2.1.1 Initial log-normal distribution 

We constructed the initial log-normal distribution by taking 
a uniform distribution (i n the log) of bu lge masses with black 
hole masses according to lHaring fc Rix (2004). Then we ap- 
plied to the bulge as well as black hole masses a log-normal 
distributed scatter with a value of ct = 0.5. Having set this 
scatter to bulge as well as to black hole masses results in a 
larger scatter in black hole and bulge masses with a value of 
~ 0.6. This scat ter is larger than the o bserv ed, present-day 
one sp ecified by ITremaine et al.l (|2002l ) and IGiiltekin et al.l 
(120091 ). We have chosen this fiducial value a posteriori as in 
our most realistic ACDM-simulation this results in a final 
observed scatter value of o" ~ 0.32 (see section [H]). How- 
ever, the scatter evolution is independent of the choice of 
the initial scatter value (see section r2.2.2[) . We start with 
an initial distribution of 580, 000 bulges with black holes. 
The black hole masses range from 2.0 < log(M./M0) < 8.0 
with a mean of (log(M,/M0)) = 5.0 and the bulges have 
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with the fitting parameters 
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Figure 2. Normalized 2D-histograni of the initial distribution 
with the same, log-normal distributed scatter {cf\^\ = 0.6) in bulge 
and black hole mass. The observed relation is shown by the black, 
solid line. The black dotted lines indicate the i-cr range of the 
initially applied scatter in our model. 
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Figure 3. Initial distribution based on the Schechter-fit for early- 
type galaxies at 2 ~ 2. An initial scatter is only applied to the 
black hole mass {cr\^\ =0.6) 



masses of 4.7 < log(MBuige/Af0) < 11.1 with a mean 
of (log(AfBuigG/-A'^0)) = 7.9. The resulting distribution of 
bulges including black holes is depicted in the 2-D histogram 
in Fig. [21 In this plot we have normalized the number of 
objects A'objccts to the maximum number of objects A'^max 
found in a black hole-bulg e mass bin. The obs erved M,- 
Mbuigc-relation according to lHaring fc Rixl l|2004h is plotted 
together with the l-cr-range of the applied scatter. 



2.1.2 Initial Schechter distribution 

Observationally, we know that the mass function o f bulges 
follows a Schechter function (e.g. I Bell et al.l 120031 ) rather 
than a log-normal distribution. Therefore we study the scat- 
ter evolution of a Schechter-shaped initial distribution of 
bulge masses. We use a fit to the measured luminosity func- 
tion (K-ban d magnitude) for red g alaxies at redshift z ~ 2 
according to ICirasuolo et al.l (|2007l ) to construct our initial 
galajcy sample, 



$(Mfc) = 0.41n(10)-'l>* • 10 



-0.4AMfc(ci-|-l) 



exp(10-°*^*^'=), 
(3) 



a 



-0.1 
-23.04 



cE)(10~^Mpc 



= 0.2 



where Mu are the absolute magnitudes in K-band. We con- 
vert the luminosity function into a mass function using the 
mass-to-ligh t ratios as function of the K-band magnitude 
according to ICappellari et al.ll2006l : 



M 



= 1., 



L VlOi«-Lfe.0 

Mfe = -2.5 log ife + 3.28 



with 



(4) 



These mass-to-light ratios were measured for a population of 
ellipticals at z = and for simplicity we assume no evolution 
with redshift. A consequence of this assumption is, that we 
obtain quite massive galaxies at 2; = 2, although the stellar 
population was not evolved completely at this time. Proba- 
bly, the mass-to-light ratio was smaller at higher z than the 
present-day value leading to smaller galaxy masses. How- 
ever, for our study it is sufficient to see the evolution of sta- 
tistical properties, which are independent of the exact choice 
of the mass-to-light ratio. We scale the resulting mass dis- 
tribution to a volume of (SOOMpc)'' and study the evolution 
of ~ 100000 bulges with black holes, considering only bulge 
masses larger than 1.6 x 10** M© (= log(A'/Buigc/M0) > 8.2). 
To keep the Schechter-distribution for the bulge masses we 
have only applied a log-normal scatter of cr = 0.6 to the 
black hole masses (see Fig. [3]). 



2.2 Depletion scenario 

For our fiducial random merging scenario (depletion case, i.e. 
without refilling the initial distribution with new galaxies) 
we use either the initial log-normal or Schechter distribu- 
tion as described above. From the initial pool we randomly 
select two objects, merge them by adding their black hole 
and bulge masses and put the merged object back into the 
pool. In the next step we again merge two objects randomly, 
but now from the new rearranged pool. This procedure is re- 
peated iteratively until, on average, every object has had one 
merger, i.e. only half of the initial objects are left over. At 
this point we define one merging generation to be completed. 
Then all remaining objects are considered as the initial pool 
for the next generation. Therefore, after the first generation 
we have N{1) = Ni-ni/'2 objects and after the n-th genera- 
tion our pool is reduced to N{n) — ^,11/2" objects. We note 
that in one generation some objects can have merged several 
times while others have not merged at all. Please note that 
if not stated otherwise, here and in the following we define 
the number of mergers by counting all mergers that occur 
for galaxies > 10 ' Mq independent of their mass ratio. 



2.2.1 Evolution of the black hole-bulge mass relation 

If we randomly merge the galaxies from the initial log- 
normal or Schechter distribution in Fig. [2] and |3l an im- 
portant consequence of the model is, as already pointed out 



p o rtant conseq r 
bv lPend ||20071 ) 



that the sample behaves according to the 
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Figure 4. Normalized 2D-iiistograins for the random merging generations 1 — 8 on tiie basis of an initial log-normal distribution, as 
shown in Fig. [2] (depletion case). The fit to the observed relation is illustrated by the black, solid line; the black dotted lines show the 
1 — cr range of the initially applied scatter. 



central-Umit-theorem (CLT). This theorem predicts that in- 
dependent of the initial distribution we always converge to- 
wards a Gaussian distribution. We see this trend already 
after only one merging generation. 

In Fig. |4]we show the evolution in the black hole-bulge 
mass plane of the log-normal distributed sample for merger 
generations n = 1 — 8. Again, the black solid line shows 
the observed, present day, M,-Mbuige-relation with the l-cr 
range of the assumed aini = 0.6 initial scatter. The relation 
is conserved during all merging generations. Note, however, 
that here we use the same initial scatter in black hole and 
bulge masses. In addition, the overall scatter decreases sig- 
nificantly with increasing merger generation. The low mass 
end of the distribution is depleted by merging whereas the 
high mass end is populated. We also see the trend that the 
scatter decreases more for more massive black holes and 
bulges than low mass systems. 

The evolution of a randomly merging initial Schechter- 
distribution of bulges is shown in Fig. [5] In this case the 
overall slope increases for small merger generations. For high 
merger generations, when the low mass end is depopulated, 
the slope again becomes similar to the initial slope and the 
relation is shifted towards larger black hole masses. The rea- 
son for the change in tilt as well as the shift is the different 
initial scatter in black hole and bulge masses. The shift to- 
wards larger black hole masses shows that having such initial 
conditions of a larger scatter in black hole masses than in 
bulge masses is quite unlikely if we want to explain the over- 
massive black holes at high redshift. However, qualitatively 
the scatter evolution is similar to the previous case. The 
scatter decreases with increasing merging generation and 
a quantitive scatter estimate is presented in the next sec- 



tion 12.2.21 We note that the initially Schechter distributed 
bulges evolve into a log-normal distribution for massive sys- 
tems as a consequence of the CLT. Still, this might not be 
a bad approximation for massive galaxies as in principle the 
Schechter distribution is a superposition of multiple Gaus- 
sians (Blanton et al. 2003). 



2.2.2 Quantifying the scatter m the black hole mass 
relation 

We characterize the scatter as the cr in a log-normal-like 
distribution: 



,f(^) = 






aV2TT 
with K = log(A/.) and n = {log{M.)). 



(5) 



Note that we use here the logarithm in base 10 'log' 
instead of the natural logarithm 'In'. This, however, only 
changes the normalization. We choose the 'log(M)' rep- 
resentation to be consistent with the observations (e.g. 
Trcmainc ct al. 2002, Giiltckin ot al. 2009). To estimate the 
scatter a for the black hole mass in the evolution of the 
M.-Mbuige-relation (Fig. |4]) we use the following method. 
For each merging generation the bulges are divided into 
different mass bins. Then for each bin we construct black 
hole mass histograms which resemble normal distributions 
which we fit with Eq. [5] to derive the scatter a. This 
method is consistent with the s catter determination in ob- 
servations (IGiiltekin et al.l 120091 ) . The fit is performed with 
a Lebenberg-Marquardt-algorithm which interpolates be- 
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Figure 5. Same as Fig.|4] but for random merging with an initial Schechter-distribution resembling 2 ~ 2 red galaxies. 



tween the Gauss-Newton algorithm and the method of gra- 
dient descent and searches iteratively for the best fit. 

Fig. [S] shows the scatter cr as a function of mean black 
hole mass (logAf.) per bin for the initial log-normal distri- 
bution (n = 0) and the eight subsequent merging genera- 
tions (n = 1...8). For the initial distribution the scatter on 
average is o — 0.6. However, it is not constant with mass 
as we have applied a log-normal distributed scatter to the 
bulge masses as well as the black hole masses. Higher merger 
generations show a decreasing average scatter for black hole 
masses larger than 10^ M0 (indicated by the dashed lines) 
as well as a decrease in scatter for larger black hole masses 
within a merging generation. This indicates a strong cor- 
relation between the scatter, the black hole mass, and the 
merging generation. 

In Fig. [7] we show the dependence of the mean number 
of mergers (A'^mcrg) on the black hole mass (log(Af.)) for 
the 8 different merging generations indicated by the colored 
lines. At the high mass end and for merger generations n ^ 7 
we find convergence towards a linear relation between black 
hole mass and mean number of mergers. The linear relation 
in Fig. [7] (red, dashed line) is given by: 



log(Armcrg) = a ■ {log(M./M0)) + b 
with a = 0.93 and b = -5.32. 



(6) 



This relation suggests, that after several merging genera- 
tions n, we can predict how many mergers a black hole of 
a certain mass must have experienced on average. E.g. a 
typical supermassive black hole of 10* Mq had about 100 
mergers, taking into account all mergers which an object 
have had during its evolution (i.e. not only mergers in the 
main branch, but all progenitors in the tree since the first 
merging generation) . Most importantly, in Fig. |8] we plot 
the scatter a versus the mean number of mergers (Afmorg) 



0.6 - 




4 5 6 7 

<log(M./Ms)> 

Figure 6. Scatter a vs. mean black hole mass (log(M,/M0)) per 
bin for different merging generations n in the random merging 
model (depletion case, initial log- normal distribution), n = is 
the initial distribution. The average values of a for black hole 
masses higher than 10^ Mq within one merging generation are 
shown by dotted lines. This shows a continuous decline in scatter 
with merger generation. 



within one bulge mass bin for the different merging genera- 
tions (indicated by diflerent colors). The decrease of scatter 
with increasing merger number is an important consequence 
of the CLT. Hence, we can derive from the CLT an analytic 
expression for the scatter decrease for an initial norm al dis- 
tribu tion as a function of the merging generation n (|Penel 
l2007f ): 



.(«) 



' CTin 



-n/2 



(7) 



where n is the generation number and aini is the initial scat- 
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Figure 7. Black hole mass (log(A/, /Af©)} as a function of the 
mean number of merger (Nmcrg) for different merging generations 
n (depletion case, initial log- normal distribution). For higher gen- 
erations n, log(A''nierg) Correlates with (log(M,/M0)). The linear 
fit is shown by the dashed, red line. 
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Figure 8. Scatter log(cr) vs. mean number of mergers 
log((Afmerg)) for different merging generations n (depletion case, 
initial log-normal distribution). The black dashed line shows the 
analytic solution according to the CLT for an initial normal 
distribution. The blue dashed line is a fit to the scatter for 
(A^mcrg) > 100, the blue one for (Nmcig) < 10. 



ter applied to black hole and bulge masses. The mean num- 
ber of mergers m of objects within one merging generation 
as a function of the generation number n can be written as 

m = 2" - 1, (8) 

Note, that for the calculation of the mean number of mergers 
we consider all merger events, which galaxies had undergone 
until this merging generation. With help of Eq. |8] we can 
rewrite Eq. [7]to get the scatter cr as a function of the mean 
number of mergers m: 

O-morg(m) ^ (Jini • (m -|- 1)~^^^. (9) 

This analytic expression is depicted in Fig. [8] by the 
black dashed line. Thereby we have made the assump- 
tion that m {— mean number of mergers per generation) 
~ (Nmerg) (= uieau uumber of mergers per bulge mass bm 
for one generation), which is a good approximation espe- 



Figure 9. Same as in Fig.|8] but based on the initial Schechter- 
distribution at 2: ~ 2. 



cially for high-mass objects. Please note, that therefore we 
will not distinguish between m and (Nmerg)- In comparison 
to merger results, the analytic solution exhibits a stronger 
decrease in scatter. This is due to the fact that t he CLT pre - 
diction assumes an initial normal distribution IjPengj 120071 ) 
whereas our sample has a log-normal distribution and there- 
fore converges towards a different scatter. The red and blue 
dashed lines in Fig. |5] are a fit to the random merging data 
assuming a fitting formula similar to Eq. |9] with the expo- 
nent a as a free parameter. 



CTin 



(m+l)-"/''. 



(10) 



For the mass range of convergence, ~ m > 100, a = 
0.61 ±0.02 (blue dashed line). The exponent a can therefore 
be used as a measure for the strength of the scatter de- 
crease. For small merger numbers (0 < m < 20) we obtain a 
weaker scatter decrease with a value of a = 0.42 ± 0.02 (red 
dashed line). This qualitatively different behaviour of the 
scatter decrease depending on the merger number range will 
be explained in section [2. 2. 3 1 If we vary the initial scatter in 
black hole and bulge mass, we obtain the same strength of 
scatter decrease within the errors, i.e. a remains unchanged. 
We have tested this for two different initial scatter values 
(Jini = 0.40 and 0.83. Even if the value to which the scatter 
converges varies, the strength of the scatter decrease is the 
same (a ~ 0.60 ± 0.02) in the limit of large m. 

In Fig.|9]we present the same scatter quantification for a 
more realistic initial Schechter distribution. Assuming con- 
vergence for ~ m > 50 or rather the small merger region 
(0 < m < 10) we get a similar value for a = 0.61 ± 0.02 or 
rather a — 0.42±0.02 (see Eg.llOp as for an initial log-normal 
distribution indicating that the strength of the scatter de- 
crease is weakly dependent on the exact choice of the initial 
distribution. Again, varying the initial scatter in black hole 
and bulge mass does not influence the value a in the expo- 
nent for large m. 



2.2.3 Difference between major and minor mergers 

According to IPena (l2007|) there is a difference for objects 
with only major or only minor mergers. He claims that major 
mergers exhibit a stronger central-limit tendency leading to 
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Figure 10. Same as Fig. [4] but galaxies had undergone only major mergers in the random merging model (depletion case). 
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Figure 11. Same as Fig. |4l but galaxies had undergone only minor mergers in the random merging model (depletion case). 



a stronger decrease of the scatter whereas minor mergers 
are mainly responsible for evolving a linear relation between 
bulge and black hole masses even if they are uncorrelated 
in the beginni ng. Going beyond the qualitative estimates 
in IPend (|2007|), we present a quantitative analysis of the 
scatter evolution for major and minor mergers. We use the 



following definitions: 

Major merger: M1/M2 < 4, Mi > M2 (11) 

Minor merger: M1/A/2 > 10, Afi > M2 (12) 

The d efinition for major mergers is consistent with IPend 
(|2007l ). In Figs.[TO]and[TT]we show the evolution of the black 
hole-bulge mass relation if we only allow major or minor 
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Figure 12. Scatter evolution for tlie depletion case as a function 
of mean number of merger for galaxies, which had either only 
major (dotted-dashed lines) or rather only minor mergers (dashed 
lines). 



galaxy mergers, respectively, for an initial log-normal distri- 
bution. In both cases the relation is conserved, however, by 
definition for minor mergers the low and intermediate mass 
range is depleted. The corresponding scatter quantification 
is shown in Fig. 1121 Indeed, there is a difference between 
majo r and minor mergers; but in contrast to the results of 
IPengj l|2007l ) in the depletion case we obtain a stronger de- 
crease of the scatter for minor mergers (a — 0.66) than for 
major mergers (a = 0.39). This result at first glance seems 
to contradict the expectations of the CLT. However, it is 
important to note, that we here consider the change of the 
scatter of black hole masses within individual bulge mass 
bins, and not over the whole population of bulge masses. 
During a merger generation a bulge mass bin has a constant 
influx and outflux of bulges due to mergers, which modifies 
the scatter behaviour with respect to the CLT. We have in- 
vestigated the scatter behaviour for major and minor merg- 
ers in the over-all distribution of bulges. We find that when 
forcing galaxies to undergo only major mergers the scatter 
even increases, whereas for minor mergers we again obtain 
a scatter decrease as expected from the CLT. This shows 
that major mergers are a very strong constraint leading to a 
violation of the CLT principle and causing the slower scat- 
ter decrease in smaller mass bins. With this behaviour we 
can deduce an explanation for the weak scatter decrease in 
the small merger number range and the stronger decrease in 
the limit of large merger number: the probability for having 
minor mergers is higher at the high mass end than for the 
low mass end, where major mergers dominate. 



2.3 Replenishment scenario 

We know from observations as well as from simulations 
that new galaxies form during the structure formation pro- 
cess. To make our simple model more realistic we now con- 
sider different replenishment scenarios. We assume again 
a certain initial number of objects A'lni with a log-normal 
or a Schechter distribution and perform the same itera- 
tive random merging procedure as described in section \2l2\ 
However, we now refill the pool after merger events with 



new objects from an external unchanged reservoir. We de- 
fine the refill-ratio Nncw/Ncvcnt to be the number of ob- 
jects added from the refill pool A'^now per number of merger 
events A'ovcnt within one merging generation. The defini- 
tion of one merging generation is the same as before but 
after the n-th merging generation the sample always con- 
tains more than N(n) = Nirii/2" objects depending on the 
refill- ratio Nnc-w / Nevsnt ■ At first we consider a refill-ratio of 
Niic-w/Ncvaut = 1, i.e. for each merger event one new object 
is added randomly from the refill pool and the total number 
of objects in the sample stays constant. We also consider 
a refill-ratio of Nnew/Ncvcnt ~ 1/3, i.e. one new object for 
every three events, motivated by ACDM simulations (see 
section [3]) . 

For the initial log-normal distribution, we consider a 
refill pool identical to the initial distribution as well as a 
pool of smaller mass galaxies with mean black hole masses 
{log(A/./MQ)) ~ 3.3 and the same initial scatter. For 
the initial Schechter distribution we use either the initial 
Schechter distribution itself as a refill-pool or we use the 
Schechter distribution containing only bulges at the low 
mass end with TTibuigo = 1.58x10® — 1.58x lO^^M©. The cases 
with smaller refill pools (lower galaxy masses) allow a more 
realistic comparison to the ACDM-simulations presented in 
sectionO In total we have four different refill-scenarios which 
will be investigated in the following: 



1. Refill-ratio 1:3 

2. Refill-ratio 1:3 

3. Refill-ratio 1:1 

4. Refill-ratio 1:1 



initial mean {Ini 1:3) 
small mean {Small 1:3) 
initial mean [Ini 1:1) 
small mean [Small 1:1) 



2.3.1 Evolution of the black hole-bulge mass relation 

A general feature of all replenishment models is that the 
scatter in the M.-MBuigc-relation is again reduced with in- 
creasing merger number. However, compared to the deple- 
tion scenario, more merger generations are needed to reduce 
the scatter by the same amount. In other words, for the same 
merger generation the scatter decrease is weaker as new ob- 
jects with a larger initial scatter are added. For small refill 
pools and large refill-ratios we find an interesting feature. 
Fig. [13] and Fig. [14] show the evolution of the relation for a 
refill-ratio of one and the small (low mass) refill pools [Small 
1:1) for an initial log-normal and a Schechter distribution. 
The contours in the plot of the first random merging gener- 
ation illustrate the distribution of the unchanged refill-pool. 
In both cases a double peak structure emerges. This is a con- 
sequence of using a low mass refill pool and a high refill-ratio. 
The low mass peak refiects the appearance of new objects 
chosen from refill-distribution whereas the high mass peak 
evolves through merging from the initial distribution, simi- 
lar to the simple case (section 12. 2p . This behaviour will be 
discussed in section [3] in more detail. 



2.3.2 Quantifying the scatter m black hole relations 

In contrast to the depletion scenario, the replenishment 
models lead to a slower decrease of the scatter in black hole 
mass. The scatter quantification for the evolution of the M,- 
-A^buigo-relation in Figs. 1 131 and 1141 (refill-ratio 1, small mean 
and initial log-normal distribution or Schechter distribution) 



10 Hirschmann et al. 




6 a 10 12 



6 8 10 12 

log{Mi,^yMo) 



6 a 10 13 

log(M,^yMo) 



6 6 10 

log(MBuiee/M©) 



Figure 13. Same as Fig. [4] but for random merging (using an initial log-normal distribution) in the replenishment scenario with a 
refill-ratio of 1:1 and a low mass refill pool. The black contours in the first merging generation indicate the distribution of the unchanged 
refill-pool. 
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Figure 14. Same as Fig. 1131 but for random merging (using an initial Schechter distribution) in the replenishment scenario with a 
refill-ratio of 1:1 and a refill pool of bulges with masses 10*Mq < A/buigo < IO^^Mq. The black contours in the first merging generation 
indicate the distribution of the refill-pool. 
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Table 1. Values of the fit parameter a for scatter decrease in 
different random merging models for a Schechter and a log-normal 
initial distribution in the limit of small (0 < m < 20) and large 
merger number (100 < rn). 

Ini. distr. Depl. Ini 1:3 Small 1:3 Ini 1:1 Small 1:1 



Schechter 












high m 


0.61 


0.56 


0.50 


0.56 


0.32 


small m 


0.42 


0.38 


0.36 


0.34 


0.54 


Log-norm 












high m 


0.61 


0.60 


0.51 


0.55 


0.31 
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0.42 


0.38 


0.34 


0.33 


0.53 
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Figure 15. Same as in Fig. [S] but for a replenishment scenario 
with a refill-ratio of 1, a refill-pool with a small mean and an 
initial log-normal distribution. The 'spike' near log(Afmerg) = 2 is 
due to the bimodality as discussed in the text. 



0.0 



i> -0.5 - 




-1.0 - 



-1.5 



-2.0 







i/a 




^ 1 




'^ 1 




11 




^^ ]- 




^ in 




n = 22 




^ ~n 





- 0.6'(m+l) 

- 0,6'(m-H) 

- 0.6*(m+l) 



log<N„ 



Figure 16. Same as in Fig. 1151 but for a replenishment scenario 
with a refill-ratio of 1, a refill-pool with a small mean and an 
initial Schechter distribution. 



is shown in Figs. [15] and 1161 The decrease of the scatter 
at low merger numbers originates from merging of new ob- 
jects from the refill-pool mainly dominated by minor mergers 
whereas the decrease of the scatter at high merger numbers 
refelects merging at the high mass end (see section 12. 2p . 
mainly dominated by major mergers. We obtain a stronger 
scatter decrease for small merger numbers (a = 0.53, red 
dashed line) than for the large merger numbers (a — 0.31, 
blue dashed line), since the probability for having minor 
mergers is higher at the small merger number end. 

We analyse the influence of refill-ratio and refill-pool 
and summarize the scatter evolution resulting from the four 
different replenishment models in table [T] for an initial log- 
normal and an initial Schechter distribtution. The typical 
errors are ±0.02. The larger the refill-ratio, the more slowly 
the scatter decreases. Keeping in mind that a larger refill- 
ratio corresponds to a larger number of new objects added 
per merger generation, we can explain this behaviour as fol- 
lows. New objects from the refill-pool have a large initial 
scatter in black hole mass, in contrast to objects after several 



Table 2. Values of the fit parameter a for scatter decrease in 
different random merging models for either only major or only 
minor mergers based on an initially log-normal distribution in 
the limit of large merger numbers (m > 100). 

Depl. Ini 1:3 Small 1:3 Ini 1:1 Small 1:1 



Major 


0.39 


0.50 


0.43 


0.54 


0.15 


Minor 


0.66 


0.65 


0.59 


0.59 


0.30 


All 


0.61 


0.60 


0.50 


0.55 


0.30 



merging generations, whose scatter has already decreased. 
Therefore, the more objects are added to the sample, the 
less the scatter decreases with merging generation. In addi- 
tion, the lower the typical mass of objects in the refill pool, 
the more slowly the scatter decreases in the limit of large 
merger numbers. That means that for large merger num- 
bers, the probability for having major mergers gets higher 
when using a low mass refill-pool. However, for low merger 
numbers the scatter decreases more rapidly with a low mass 
refill-pool. In this range, minor mergers become more likely 
because of the low-mass refill-pool, leading to a stronger 
scatter decrease than in the major merger dominated re- 
gion (i.e. convergence region of large merger numbers). Al- 
together, a small-mass refill-pool leads to higher probability 
of having minor mergers in the small merger number re- 
gion than at the limit of large merger numbers, where ma- 
jor mergers dominate. Furthermore, as shown in table[TJ the 
choice of the initial distribution does not change the scatter 
decrease a (for small as well as large m). We expect this from 
the CLT as the distribution quickly evolves into a normal 
distribution no matter which initial distribution is used. 



2.3.3 Difference between major and minor mergers 

Since IPengj (|2007l ) considered in his study a replenishment 
scenario with a refill ratio of A'ncw/A'^ovont = 1 using the ini- 
tial distribution as the refill pool, we will also investigate 
the difference between major and minor mergers for differ- 
ent replenishment models. The fitted slopes for the scatter 
evolution for major and minor mergers are summarized in 
table E33] with an error of ~ ±0.02. Note that for the case 
refill-ratio 1 : 1 with a low-mass refill-pool, we obtain a 
larger error of ~ ±0.05, since here the black hole mass his- 
tograms are not fitted well by a Gaussian function anymore. 
As in the depletion case (see sect ion [2. 2. 3 1) we see a stronger 
scatter decrease for galaxies undergoing only minor mergers 
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than for major mergers. We argue as above that this devia- 
tion from the CLT is most likely due to dividing the bulge 
masses in different mass bins which suffer a constant influx 
and outflux from bulges during each merger generation. And 
furthermore, forcing galaxies to undergo only major mergers 
is a very strong constraint leading to a violation of the CLT. 



3 COMPARISON TO MERGING IN 
ACDM-SIMULATIONS 

So far, the influence of different idealised random merging 
models on the evolution of the M,-Mbuigc-relation and the 
corresponding scatter in black hole mass has been discussed. 
In this section we investigate a more complex and astrophys- 
ically motivated model based on merger trees from dark mat- 
ter simulations following structure formation in a ACDM 
universe. 



3.1 Simulation setup and merger tree 
construction 

We have simulated a comoving periodic box with L = 
100 M pc box length and 5 12^ particles using GADGET2 
code (jSpringel et al.l [200l|). For this simulation we as- 
sume a ACDM cosm ology based in the WMAP3 (see e.g. 
ISpergel et al.l ((2003)) measurements with ag, — 0.77, Q,m = 
0.26 ^A = 0.74 and h = //o/(100 kms'^) = 0.72 (see 
also iMoster et al.l (120091 ) for a first analysis of this sim- 
ulation). The simulation was started at z = 43 and run 
until z = G with a fixed comoving softening length of 
2.52 /i^^kpc. Starting at an expansion factor of a = 0.06 
we have halo catalogues for 94 snapshots until z = Q sepa- 
rated by Aa = 0.01 in time. The mass of one dark matter 
particle is 2 x l{f Mq / h. To identify dark matter haloes at 
every snapshot we use a FOF algorithm with a linking length 
of 6 = 0.2. In a second step we extract sub haloes of every 
FOF halo using the SUBFIND algorithm (|Springel et all 
[200l|). This halofinder identifies overdense regions and re- 
moves gravitationally unbound particles. This way we split 
the FOF group into a main or host halo and its satellite 
halos. The sizes and virial masses of the main halos (i.e. 
the most massive SUBFIND halos) are determined with a 
spherical overdensity criterion. The miminum halo mass is 
set to 20 particles (4 x lif Mq/K). To create the merger tree 
for the main halos (satellite halos are not considered) we 
co nnect halos betwe e n the snapshots as described in detail 
in iMaulbetsch et al.l (|2007l ) with a few modifications which 
we will explain in the following. The branches of the trees 
for 2 = halos are constructed by connecting the halos to 
their most massive progenitors (MMP) at previous snap- 
shots. Thereby halo j with Uj particles at redshift Zj with 
the majcimum probability p{i,j) is chosen to be a MMP of 
halo i containing Ui particles at redshift Zi (where j < i). 
The probablity p{i,j) is defined as 



P(hj) = 7^^ With 

Uov = riiizi) n nj{zj) and 
nm.ax{i,j) = max{ni{zi),nj{zj)) 



(13) 



Here, nov is the number of particles found in both halos 
and Umax is the particle number of the larger halo. We re- 
move 'fake' haloes which exist only within one timestep and 
have no connection to any branch. The low redshift ends 
of the branches are then checked for mergers. A halo j is 
assumed to merge with halo i, if at least 50% of the parti- 
cles of halo j are found in halo i. In case of a merger the 
branches are connec ted. For a pr oper connection we apply 
the 'split '-algorithm IjGenel et al.ll2008 l to prevent double or 
multiple counting of merger events within a tree. The 'split' 
algorithm was shown to produce more reliable merger rates 
which is also important for our study. 



3.2 Evolution of the relation between black hole 
mass and galaxy properties 

We consider two different possibilities to populate dark mat- 
ter haloes with black holes and bulges. We either apply black 
hole masses or bulge masses directly to dark matter halos 
and calculate the corresponding miss ing quantity u s ing th e 
the black hole-bulge mass relation in lHaring fc Rixl (|2004f ). 



3.2.1 Using a M, -Md m -relation 

iFerraresd (12003) proposed a relation between the mass of 
the central black hole M. and the mass of the dark matter 
halo AIoM of the form 



M. 
108 Mp 



0.1 



Mb 



1012 Af( 



(14) 



We want to point out that only the M,-Mbuigc-relation is 
directly observable and therefore assumed to be more fun- 
damental than the M,-MDM-relation. However, the M,- 
AfoM-relation h a s bee n supported by e.g. recent results of 
IBooth fc Schavd (120091), where the authors obtain a similar 
relation between black hole and dark matter halo mass using 
numerical simulations (GADGET HI) with self-consistent 
black hole growth, which are tuned to match the relations 
between black hole mass and galaxy s tellar properti e s. Fur - 
thermore there are observations fro m ICroom et al.l (|2007|) ; 
IYu fc Lul (|2008l ): IWhite et al.l (|2008h which confirm the exis- 
tence of a relation between black hole and dark matter halo 
mass using observations. However, this relation is only valid 
for z — 0. Since we seed the black holes at higher redshift, we 
modify expression Eq. [T?] in order to maintain the relation 
at 2 = 0. We ass ume the same derivation as it is described in 
IFerraresd (|2002r ). Therefore we first use the virial velocities 
Vvir from our simulations calculated by SUBFIND. This way 
we obtain v^ir as a function of Mdm and redshift z. Assum- 
ing that v-uir ~ Vc we use then the relation between circular 
velocity Vc and velocity dispersion cjc and the one between 
velocity dispersion cjc and black hole mass M, . So we can 
rewrite equation 1141 



M. = 3.12 X 



Vc{M, 



DM, Z) 



1.19 \ 4.58 



200 km/s 



10'' Mq 



(15) 



Note, that the relation between circular velocity Vc and ve- 
locity dispersion CTc is derived from observations of spiral 
galaxies at 2; = 0. Presumably at higher redshifts dissipa- 
tion of baryons has a higher influence than for 2 = 0. This 
could lead to larger velocity dispersions and therefore also 
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Figure 17. Evolution of the M,-Mbuige-rclation through merging only in ACDM-simulations using A/.-AfoM-rclation. 



to larger black hole masses. According to relation [15] we can 
populate halos with masses extracted from the dark matter 
simulations with central supermassive black holes. As we 
are mainly interested in massive galaxies where gas physics 
is assumed to be less important (at least at low redshifts e.g. 
Dekel fc Birnboimll2008( : iNaab et al.ll2006l : iKhochfar fc Silkl 
2009|) we only investigate merger trees starting at z = 3 for 
halos more massive than lO^'^M© at 2 = 0. At the high red- 
shift end of every branch starting at z = 3 we populate the 
most massive progenitors of the selected 2 = halos with 
black holes according to Eq. 1151 Additionally we add a log- 
normal distributed scatter to the black hole masses with the 
fini = 0.6 (see section [2}. For the subsequent growth of the 
black holes we only take dark matter halo mergers and the 
corresponding merger of their black holes into account. 

To allow a comparison with the random merging cases 
shown previously, the evolution of the M.-Mbuigc-relation 
is shown in Fig. 1171 We have chosen black holes in dark 
matter haloes according to Eq. [15] and then calculated cor- 
res ponding bulge mass es by taking the median relation 
of JHaring fc Rixl (|200J). Once the dark halos have been 
populated with these bulge masses we apply a scatter of 
sigma=0.6 to the black hole masses for each bulge mass. 
We again follow the growth process only via merger events. 
The high mass end of the relation is shifted towards larger 
black hole masses as we have only applied an initial scatter 
to black hole masses but not to the bulge masses. Similar 
to the simple models investigated before again the scatter 
decreases with time. We also see a double peak structure 
at low redshifts z < 0.4 similar to the replenishment ran- 
dom merging model with a low mass refill pool (section l2.3p . 
This is a consequence of the conditional mass function. We 
only investigate halos under the condition that they merged 



into 2 = halos with masses larger than IO^^Mq (see e.g. 
ISomerville et al.ll2000l ). 



3.2.2 Using a galaxy population model 

Alternatively, to populate dark matter halos with black 
holes, we can use a fitting function that relates host dark 
halo masses to stellar masses of galaxies to assign to ev- 
ery dark matter halo mass a galaxy stellar mass. There ex- 
ist many studies which li nk the distribution of galaxies to 
that of dark matter halos (^van don Bosch et al. 2003, 200^ 
jMandclbaum ct al. 2006; Mostcr ct al. 2009; Guo & Whit3 
|2009|) . Here we take the fitting formula from JMoster et al.l 
(J2003) . They assume that every host halo contains exactly 
one central galaxy and - as a constraint from the observed 
galaxy mass function - that the stellar mass to dark mat- 
ter halo mass ratio Mt/MoM first increases with increasing 
mass, reaches a max imum and then decreases again. Hence 
JMoster et al.l (|2009l ) ado pt the following p arametrization, 
similar to the one used in I Yang et al.l l|2003l '): 



M4Mc 



Me 



M, 
Mdm 



(MdmY 

\ Ml ; 



f Mdm \ 
\ Mr ) 



(16) 
Basically, this parametrization is set to reproduce many ob- 
servations, as the galaxy mass function or clustering. Choos- 
ing a redshift parametrization for each of the parameters in 
Eq. 1161 they can predict the galaxy to dark matter mass 
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Figure 19. Scatter in black iiole mass versus mean black hole 
mass (log{A'/, /Mq)) at different redshifts (different colors) in the 
ACDM-simulation assuming black hole seeding according to the 
M.-M^i^Vf-relation. The horizontal lines indicate the average scat- 
ter for black holes more massive than ^ 10^ Mq. The scatter for 
massive black holes continously decreases towards lower redshifts. 



Figure 20. Relation between mean black hole mass 
(log(A/, /A/q)) and mean number of mergers log(A''niorg) at 
different redshifts (different colors) in the ACDM-simulation 
assuming black hole seeding according to the M.-M^ij^f-relation. 
The best fit at z=0.05 is indicated by the red dashed line. For 
comparison the dashed black line indicates the best fit for the 
random merging depletion model (Fig.0. 



ratio at any redshift: 

log Ml (2) =logMo-(2-l-l)'' 
M. \ , , / M. 



Mr 



7(2)=7o-(2 + ir 



(17) 

(18) 

(19) 
(20) 



with logA/o = 11.88, ^ = 0.019, {M*/Mdm)z=o = 0.0282, 
V = -0.72, 70 = 0.556, 71 = -0.26, /3o = 1.06 and /?i = 
0.17. 

To populate the galaxies with black holes, we assume 
for simplicity that all stars are in the spheroidal component 
of the galaxy (M, » Mi,uige)- Using the Af,-M(,uige-relation 
(JHaring fc Rixll2004l ) we apply to each galaxy mass a black 
hole mass with the same initial scatter as already used be- 
fore (a = 0.6). If we then take into account the growth 
of black holes and galaxies through merging according to 
our ACDM-simulations, we obtain an evolution of the M,- 
Afbuige-relation as shown in Fig. 1181 Again, we see the same 
effect as already described in subsubsection [3T2 . 1 1 a decreas- 
ing scatter with time together with an evolving double peak 
structure. 



3.3 Quantifying the scatter in the black hole mass 
relation 

For the evolution of the M.-Mbuigc-relation shown in Fig. 
1171 (the model based on the relation between dark halo mass 
and and black hole mass) we performed the same scatter 
quantification as before for the different random merging 
models. Fig. [!§] shows the evolution of the scatter in black 
hole mass. At the high mass end, which we are focusing 
on here, the scatter decreases towards lower redshift. Fur- 
thermore - similarly to random merging - in Fig. [20] we find 
convergence towards a linear relation between the logarithm 
of the mean number of mergers and the logarithm of black 
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Figure 21. Scatter log(cr) as a function of mean number of merg- 
ers log(A'^morg) for different time steps in the ACDM-simulation 
assuming black hole seeding according to the M.-M^i jvf-relation. 
The blue dashed line shows the fit to the best-matching random 
merging model (refill-ratio 1/3 and low-mass refill-pool), the black 
dashed line corresponds to the analytic case. 



hole mass at low redshift. We fit this relation according to 

log{iV„,erg) = a ■ {log(M./M0)) + b (21) 



with a = 0.46 and b - 



-1.39. 



shown by the red dashed line in Fig. 1201 For comparison the 
dashed black line illustrates the random merging case (de- 
pletion). As expected lower mass black holes (M. < 10* M©) 
experience more mergers than predicted for the simple de- 
pletion case, whereas at the high mass end galaxies undergo 
less mergers than in the depletion model. 

Finally the scatter a is plotted as a function of the 
mean number of mergers (A^'mcrg) in Fig. 1211 We find the 
same qualitative behaviour as for all random merging mod- 
els. The scatter decreases with increasing number of merger 
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Figure 18. Evolu tion of the M.-Mbuige-relation through merging only in ACDM-siniuIations using a 
l|Moster et al.ll2009h . 



;alaxy population model 



events. However, the scatter decrease in ACDM-merging is 
weaker. The fit parameter a — 0.30 ± 0.03 is smaller than 
for most of the random merging models at the limit of large 
merger numbers. Note that for the case where we assign 
black holes to dark matter halos using the galaxy popu- 
lation model, as described in subsection 13.2.21 we get the 
same qualitative scatter behaviour. The only difference is 
that because there are fewer objects at the high mass end, 
the resulting scatter value is connected with a larger error 
than for the black hole population described in subsection 
13.2.11 In comparison to the best-matching random merg- 
ing scenario - a refill-ratio of 1/3 and a low-mass refill-pool 
- ACDM-merging leads to a clear difference in the scatter 
decrease (a = 0.51 ± 0.02). However, in the case of ACDM- 
merging we have only a maximum number of mergers of 
~ 60. Therefore, the scatter decrease in CDM-merging is 
consistent with the best-matching random merging model 
only in the low merger range (a — 0.34 ± 0.02), indicated 
by the blue, dashed line. This shows that the scatter evolu- 
tion in ACDM-merging can be well approximated - quanti- 
tatively and qualitatively - by the simple model of random 
merging without any structure formation model. 



3.4 Evolution of the black hole mass function 

In this section we investigate the evolution of the black hole 
meiss function in our merger-driven model to test, if merging 
only provides an adequate description for black hole growth. 
Fig. 1221 shows the black hole mass function for different red- 
shifts (different colors) assuming an initial log-normal scat- 
ter in black hole mass at z = 3 and a growth of black 
holes via mergers only. The solid lines indicate the evolu- 
tion of the black hole mass function assuming seeding ac- 



cording to the M,-MDM-relation (JFerraresel l2002f ) and the 
dashed lines show the black hole mass function b ased on the 
galaxy population model of lMoster et al.l (|2009l ). In the left 
panel we show the local observed black hole mass function 
(black triangles) derived from correlations between black 
hole mass and bulge lum inosity or stellar velocity disper- 
sion ([Marconi et al.ll2004l '). This is in disagreement with the 
2 = model prediction (red line), since we underestimate 
black hole masses. This implies that growth by gas accretion 
might be an important contribution to the overall formation 
process of black holes and should not be neglected. Further- 
more, if we seed the black holes only at z = 1, the disagree- 
ment with observations is - as expected - less pronounced 
than for z = 3-seeding. But again, the black holes are not 
massive enough to reproduce the observations. This shows 
that even from z = 1 till z = merging only seems to be an 
unsufBcient description for black hole growth. However, we 
want to point out most importantly that in this study we are 
not trying to fit the present-day black hole mass function. In 
fact we want to show how the scatter evolution in black hole 
mass would be affected by merger events only. The possible 
influence of additional gas physics will be discussed in the 
next section. 

The right panel of Fig. [22] shows the black hole mass 
function for the two different seeding mechanisms only at 
z = 0. However, here we have considered only galaxies that 
have undergone a merger with a merger ratio smaller than 
1 : 10 in the last 100 • 10'^ yrs. Therefore, they can be as- 
sumed to be in the active phase during this time. Our results 
deviate from observed values of the present-day black hole 
mass function fo r a sample of active galaxies (8500 objects 
from SDSS DR4. lGreene fc Holl2007h . On the one hand, this 
could be a consequence of underpredicting the overall black 
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hole mass function as shown by the left panel, assuming 
the correct number of merger events. Alternatively, it might 
imply that merger events are not frequent enough to give 
an explanation for observed active galaxies and accretion is 
needed even without any merger event. However, if we as- 
sume, that also lower mass ratios (< 1 : 100) can trigger 
the activity of galaxies, we obtain a good agreement with 
observations. Also raising the time for one duty cycle would 
lead to a better consistence with observations. This could 
be justified by the fact, that within our assumed time for 
one duty cycle there could exist more active objects if we 
would consider galaxy mergers, which happen later in time 
than mergers of isolated dark matter halos (which are taken 
into account in this study). 



4 CONCLUSIONS AND DISCUSSION 

In this paper we investigated the evolution of the intrinsic 
scatter in black hole mass under the assumption that the 
black holes only grow by mergers with other black holes 
during galaxy mergers. For many different idealized random 
merging models with (replenishment) and without (deple- 
tion) refilling from an external galaxy pool we find the fol- 
lowing general results: 

(i) The evolution of the black hole distribution can be well 
described within the framework of the central limit theo- 
rem (CLT). Independent of the initial distributions, e.g. log- 
normal or Schechter, of black holes the distribution always 
evolves into a normal distribution after a few merger gen- 
erations. We consider all mergers independent of the mass 
ratio that galaxies > 10*'^ M© experience. 

(ii) For all random merging models we found a decreas- 
ing scatter a with increasing merging generation n and with 
increasing merger number m. As a consequence of the mass 
built-up during the merger events the scatter also decreases 
with increasing black hole mass. Motivated by the CLT, we 
can approximate the scatter dependence on the mean num- 
ber of mergers by 

Crmcrg(m) « (Tini ■ (m -f 1)"°'^. (22) 

Here the exponent a is a measure of the strength of the 
scatter decerease. For the different random merging models 
we find 0.30 < a < 0.61 for a large number of mergers (m > 
100) independent of the initial scatter applied to black hole 
and bulge masses. In general, replenishment models show a 
weaker scatter decrease. 

(iii) Considering either only major or only minor mergers 
for galaxy growth we found that minor mergers lead to a 
much stronger scatter decrease than major mergers; hence 
the smaller the mass ratio of merger events, the more r apidly 
the sc atter decreases. This is in contrast to findings of lPeng 
l|2007t l. 

(iv) For different replenishment models we found that the 
higher the refill-ratio and the smaller the typical mass of 
black holes in the refill-pool, the more slowly the scatter 
decreases. 

Studying the effect of merging according to current 
structure formation models in ACDM-simulations, we find a 
qualitatively similar behaviour. The scatter decreases with 
the number of mergers (a = 0.3), and as a consequence it 



also decreases with cosmic time. This is also quantitatively 
consistent with the best-matching random merging model 
(refill-ratio 1/3 and low mass refill-pool) at least for the 
limit of the low merger number range (a = 0.34), since in 
ACDM-merging the most massive galaxies experience only 
~ 50 — 60 merger events. Therefore, the scatter evolution 
in ACDM-merging can be well approximated by a simple 
model assuming random merging of galaxies. 

From the above results we can draw some implications 
about recent observations of high redshift black holes: 

• For a simple merger driven growth of black holes we 
predict that the scatter in black hole mass must have been 
larger at higher redshift. Assuming an initial scatter of 
0.6 at z = 3 the ove r -massive black ho les investiga ted by 
[Schramm et~ai] (|2008l ). |Peng et all l|2006l ) and McLure et all 

(|2006r ) would be within the 2 — a range of this large initial 
scatter. If these objects had on average 50 — 60 dry merger 
events then we also obtain the present-day scatter v alue of 
~ 0.31 for massive ellipticals (|Giiltekin et al.ll2009l ). This 
shows that the observations of over-massive black holes at 
high redshifts are consistent with a population of galaxies 
that has a large scatter in black hole mass (~ 0.6) at high 
redshifts being reduced through subsequent merging. 

• A further important advantage of our model is that 
we are not only able to explain over-massive black holes 
at high redshifts, but we can also ac count for the observed 
undcr-massivc black holes at 2; = 2 ^Alexander et al.l 120081 : 
[Shapiro ct al. 2009). These objects would be lying within 
the 2 — (T range of the scatter a± z = 2 [a ~ 0.5) assuming 
an initial scatter of 0.6 at z = 3. However, we want to point 
out, that models assuming an evolution of the M. — Mbuige 
relation (from higher towards lower black hole masses with 
decreasing redshift) can be an explanation for over-massive 
black holes at high redshifts, but not for the observed under- 
massive ones. 

• We can explain th e scatter evolution o f model predic- 
tions from GALFORM (|Malbon et aLl bOO?'): the decreasing 
scatter with increasing black hole mass (see Fig. [l|. This is 
the consequence of the merger driven growth of black holes. 
More massive objects have undergone more merger events 
than less massive ones; this higher merger number leads to 
a stonger decrease in scatter. However, observations show 
a larger s catter in black hole masses, especially at the high 
mass end JTremaine et al.ll2002l . ICiiltekin et al.ll2 009'). That 
might indicate that in the model of iMalbon et al.l (|2007l ) ei- 
ther the number of merger events for massive objects was 
overestimated (unlikely) or that they started with a too 
small scatter in black hole mass at high redshift. 

• Our finding of a decreasing scatter with increasing 
black hole mass as a consequence of subsequent merging 
are also consistent with the results of[Gultckin ct al. (200!|), 
who distinguish between two subsamples of ellipticals and 
non-ellipticals. For ellipticals they find a smaller scatter in 
black hole masses than for non-ellipticals. Referring elliptical 
galaxies mainly to the high-mass end and non-elliptical ones 
mainly to low-mass end, these observations are also consis- 
tent with our results. Furthermore, ellipticals presumably 
had more mergers. 

The aim of this paper is to investigate the scatter evo- 
lution of supermassive black hole for dry merger driven 
growth. Despite the interesting insights we have neglected 
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Figure 22. Left; Evolution of the black hole mass function through ACDM-merging using using Af,-A/i5j(/-relation (solid lines) a nd using 
the ga laxy population model (dashed lines). The triangles correspond to the local black hole mass function according to Marc oni et alj 
||200J). Right: Same as in the left panel, but only for z = 0. However, considered are only objects, that had undergone a merger event 
within the last 100 ■ 10^ yr s with a merger rat io smaller than 1 : 10. The diamonds show the observed, present-day black hole mass 
function for active galaxies JGreene &c Hdl2007r i. 



another important physical growth mechanism, e.g. the ac- 
cretion of gas onto black holes. The importance of an addi- 
tional growth mechanism is also confirmed by the compar- 
ison of the black hole mass function resulting from merg- 
ing only with observations at z — 0: we under-estimate the 
masses of black holes. Furthermore, if we consider only ac- 
tive galaxies that have undergone a merger event in the last 
100-10® yrs with a merger ratio smaller than 1 : 10, we again 
deviate from observations, which could be a consequence of 
underpredicting the over-all black hole mass function. Al- 
ternatively, it might indicate that merger events are not fre- 
quent enough to be an explanation for the observed active 
galaxies. This implies that there has to be additional gas 
accretion even if no merger event happened. 



F urthermore, according to the Soltan argument (|Soltanl 
Il982l ) accretion on to massive black holes is the dominant 
source of energy produced by quasars (L ~ M,c^), which 
presumably indicates that the mass in black holes is mainly 
generated by accretion of gas. This gas accretion might be 
triggered by galaxy mergers and their ability to drive large 
amounts of gas to t he centers of the galaxies an d possi- 
bly feed a black hole ( Springel et al. 2005; Di Mat teo et al.l 
l2005l : [Johansson et all l2009l :" lHopkins et all |2008| ). It was 
shown in hydrodynamical galaxy-galaxy merger simulations 
with self-regulated black hole growth that after the active 
phase of the quasar, star formation and black hole growth 
are quenched by gas heating through energy release by 
the active quasar. Following the merger, the quasar host 
becomes a red and dead, massive elliptical galaxy with- 
out any further black hole growth through gas accretion 
l| Johansson et al.ll2009l ). For these galaxies the further evo- 
lution of their black holes is mainly dominated by dry 
merger events. However, gas-rich mergers and cold accre- 
ation flows seem to be the dominant growth mode for mas- 
sive high redshift black holes, as th e galaxies are more 
gas rich (JKhochfar fc Silkll2006l . |2009| ). The quasar activ- 



ity p eaks at z ^ 2 ~ 3 (e.g. iHasinger et al.l l2005l : lUedal 
[200a). Therefore, this phase is most likely to be respon- 
sible for the scatter in black holes at this redshift (corre- 
sponding to th e initia l scatter in our models) . According to 
IHasinger et al.l (|2005l ) the emissivity of the high luminos- 
ity objects (log(L/L0) ^ 45) drops by ~ 10^ ergs"^Mpc"^ 
between z ~ 2 and 2 ~ 0. Therefore, the subsequent evolu- 
tion, at least for these massive black holes, might be driven 
mainly by dry mergers as investigated in this paper. In 
addition, there is theoretical evidence that for z ^ 2 gas 
is more efficiently heated by accretion shocks associated 
with gravitational heating in massive ( Mdm ^ 10"^^ Mq) 
galaxy halos (Khochfar & Ostr ikcr 2008; Dekel fc BirnboimI 
I2OO8I : iNaab fc Ostrikeij l2009l : iBirnboim et al.l |2007| ). sup- 
pressing further gas cooling and star formation as well 
as accretion onto the central supermassive black hole. 
Thus, 'red and dead' massive spheroids evolve starting 
at z ~ 1, where again merging may dominate the fur- 
ther growth of these b lack holes IjDekel fc BirnboimI 120061 : 
iKhochfar fc Silkl 2009). From disk galaxy merger simula- 
tions in I Johansson et al.l l|2009l ). we know that overmassive 
black holes lying above the black hole-bulge relation do not 
evolve onto the relation considering only gas accretion driven 
growth. A solution for this problem could be that the evo- 
lution of overmassive black holes onto the relation might be 
caused mainly by growth through merger events leading to 
a scatter decrease in black hole mass as shown in this paper. 
That would also indicate, that gas accretion does not play 
the most important role in the evolution process of overmas- 
sive black holes, at least at late stages, where dry merging 
is the dominating process. 

Intentionally, in this paper we kept the physics very 
simple and focused only on the growth through merging in 
order to understand the influence of merging on the evolu- 
tion of the scatter in black hole mass. However, gas accretion 
is - even for z < 2 - an important growth channel. The im- 
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pact of gas accretion on the evolution of the scatter in the 
M, — Mbuige relation is complicated and so far not easy to 
assess. T o get an idea for such an influence, we refer the 
reader to [Johansson et alj (|2009|). Using merger simulations 
of disk galaxies with gas (3 • 10® M© < Mbuige < lO^M©) 
they show two possible scenarios: Starting with galaxies on 
the relation leads to an increase of the scatter which gets 
even larger for higher gas fractions. However, considering 
initially galaxies below and above the relation causes a de- 
crease of the scatter. The decrease seems to get stronger for 
lower gas fractions. In order to really understand the effect 
of gas physics we definitely need further investigations of 
this process in a statistical sense. 
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